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ABSTRACT 

We revisit the atmospheric properties of the extremely hot Jupiter WASP-12b in hght of several 
new developments. First, new narrowband (2.315 /xm) secondary eclipse photometry that we present 
here, which exhibits a planet/star flux ratio of 0.45 % ± 0.06 %, corresponding to a brightness temper- 
ature of 3640 K ± 230 K; second, recent Spitzer/IRAC and Hubble/WFC3 observations; and third, 
a recently observed star only 1" from WASP-12, which has diluted previous observations and which 
we further characterize here. We correct past WASP-12b eclipse measurements for the presence of 
this object, and we revisit the interpretation of WASP-12b's dilution-corrected emission spectrum. 
The resulting planetary emission spectrum is well-approximated by a blackbody, and consequently 
our primary conclusion is that the planet's infrared photosphere is nearly isothermal. Thus secondary 
eclipse spectroscopy is relatively ill-suited to constrain WASP-12b's atmospheric abundances, and 
transmission spectroscopy may be necessary to achieve this goal. 

Subject headings: infrared: stars — planetary systems — stars: individual (WASP-12, Bergfors-6) — 
stars: multiple — techniques: photometric — techniques: spectroscopic — eclipses 



1. INTRODUCTION 

Transiting extrasolar planets allow the exciting possi- 
bility of studying the intrinsic physical properties of these 
planets. The latest new frontier to emerge is the detailed 
study of molecular chemistry in the atmospheres of these 
planets, many of which exist in intensely irradiated en- 
vironments. Recent years have seen rapid strides in this 
direction, with measurements of precise masses and radii, 
detection of secondary eclipses and phase curves and 
the start of ground-based spectroscopy (Redfield et al. 
2008; Swain et al. 2010; Bean et al. 2010). Based on 
observed day/night temperature contrasts (e.g.. Cowan 
& Agol 2011), atmospheric circulation patterns (Knut- 
son et al. 2009), and atmospheric chemistry (Stevenson 
et al. 2010; Madhusudhan et al. 2011) these planets' at- 
mospheres are likely to be quite different from anything 
previously known. 

1.1. Introducing the WASP-12 System 

A prime example is the transiting Hot Jupiter WASP- 
12b, which is one of the largest and hottest transiting 
planets known (Hebb et al. 2009; Chan ct al. 2011; Ma- 
ciejewski et al. 2011). The planet is significantly over- 
inflated compared to standard interior models (Fort- 
ney et al. 2007), though its radius and age can be ex- 
plained by an appropriate dynamical history involving 
an initially eccentric orbit and subsequent interior dis- 
sipation of tidal torques (Ibgui et al. 2011). Radial ve- 
locity measurements associated with the initial transit 
discovery and the first occultation observation both sug- 
gested WASP- 12b had a nonzero eccentricity (Hebb et al. 

^ Max-Planck Institiit fiir Astronomie, Konigstuhl 17, D- 
69117, Heidelberg, Germany; ianc@mpia.de 

Department of Physics & Astronomy, University of Califor- 
nia Los Angeles, Los Angeles, CA 90095, USA 

3 Lowell Observatory, 1400 West Mars Hill Road, Flagstaff, 
AZ 86001, USA 

* Subaru Telescope, National Astronomical Observatory of 
Japan, 650 North A'ohoku Place, Hilo, HI 96720, USA 



2009; Lopez-Morales et al. 2010). However, subsequent 
orbital characterization via timing of secondary eclipses 
(CroU ct al. 2011; Campo ct al. 2011; Cowan et al. 2012) 
and further radial velocity measurements (Husnoo et al. 

2011) set an upper limit on the eccentricity of ^ 0.03(lcr). 
The 2.315 /zm narrow band eclipse we present here is also 
consistent with a circular orbit. 

Due to its close proximity to its host star the planet is 
thought to be significantly distorted and may even be un- 
dergoing Roche lobe overfiow (Li et al. 2010). Such over- 
flow, if verified, would be the first evidence of the tidal 
inflation instability (Gu et al. 2003). Possible evidence 
for the overflow scenario has come (1) from HST/COS 
UV spectra taken during transit (Fossati et al. 2010), 
which show tentative evidence of a deeper transit with 
earlier ingress than observed in the optical (Hebb et al. 
2009), (2) from a tentative detection of an extended Ks 
band secondary eclipse duration (Croll et al. 2011), which 
could be interpreted as an opaque accretion stream or 
disk, and (3) from Spitzer/IRAC phase curve observa- 
tions of WASP-12b, which detect ellipsoidal variations 
from the planet at 4.5 /im at a level consistent with a 
planet filling (or overfilling) its Roche lobe (Cowan et al. 

2012) . However: (1) there is no evidence for an ex- 
tended occultation duration in Spitzer/IRAC observa- 
tions (Campo ct al. 2011) or in the 2.315 /im narrowband 
eclipse we present here; (2) degeneracies between ellip- 
soidal variations, thermal phase variations, and instru- 
mental systematics prevent an unambiguous determina- 
tion of WASP-12b's geometry from the Spitzer observa- 
tions (Cowan et al. 2012); and (3) recent HST/WFC3 
secondary eclipse spectroscopy suggest that WASP- 12b 
is not substantially distorted (Swain ct al. 2012). 

WASP-12b is intensely irradiated by its host star, mak- 
ing the planet one of the hottest known and giving it 
a favorable {> 10~^) NIR planet/star flux contrast ra- 
tio; its atmosphere has quickly become one of the best- 
studied outside the Solar System. The planet's large size. 
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low density, and high temperature motivated a flurry of 
optical, (Lopez-Morales et al. 2010), NIR (CroU ct al. 
2011), and mid-infrared (Campo ct al. 2011) secondary 
eclipse photometry has been interpreted to reflect an at- 
mosphere with an unusual carbon to oxygen (C/0) ratio 
greater than one (Madhusudhan et al. 2011). Subsequent 
ground-based observations (Zhao et al. 2012; Crossfield 
et al. 2012) and the recent WFC3 1.1-1.7 ^m spectrum 
(Swain et al. 2012) are consistent with these earlier mea- 
surements and the C/0> 1 model, but the 2.315 /im 
eclipse we present here is inconsistent (at > 3a) with 
such models. In addition, under the so-called "null hy- 
pothesis" (i.e., a spherical planet) of Cowan et al. (2012) 
the IRAC 4.5 /itm secondary eclipse is significantly deeper 
than the previous measurement (Campo ct al. 2011), sug- 
gesting less absorption by CO and weakening the case for 
a high C/0 ratio. 

As yet transmission spectroscopy (which determines 
atmospheric opacity at a planet's limb via multi- 
wavelength transit measurements Seager & Sasselov 
2000) has so far been limited for this system. Opti- 
cal transit measurements show some disagreement (Hebb 
et al. 2009; Chan et al. 2011; Macicjcwski ct al. 2011), 
which makes interpretation difficult. Spitzer/IRAC tran- 
sit observations suggest that WASP-12b's radius may be 
greater at 3.6 /xm than at 4.5 fim (Cowan et al. 2012), but 
only if the planet is much more prolate { Riong / Rp = 1-8) 
than suggested by WFC3 observations (3cr upper limit 
of 1.7; Swain ct al. 2012). Under the null hypothesis of 
Cowan et al. (2012), the transit radius is larger at 4.5 /im 
(as expected from atmospheric models). WASP-12b's 
low density and high temperature ensure that this planet 
will continue to be a target for future efforts in this direc- 
tion; if (as we suggest) the planet's atmosphere is in fact 
nearly isothermal at the pressures probed in secondary 
eclipse, transmission spectroscopy may be the only hope 
for constraining WASP-12b's atmospheric composition. 

Thus significant uncertainties remain in the interpre- 
tation of the current ensemble of atmospheric measure- 
ments. At the moment this situation is typical even for 
the best-characterized systems (Madhusudhan & Scagcr 
2010) because (a) broadband photometry averages over 
features caused by separate opacity sources and (b) at- 
mospheric models have many more free parameters than 
there are observational constraints. When properly cali- 
brated, spectrally resolved measurements can break some 
of these degeneracies. Such results can test the interpre- 
tation of photometric observations at higher resolution, 
and can more precisely refine estimates of atmospheric 
abundances, constrain planetary temperature structures, 
and provide deeper insight into high-temperature exo- 
planetary atmospheres. These goals provided the mo- 
tivation for our earlier ground-based spectroscopy of 
WASP-12b (Crossfield ct al. 2011) and serve as the im- 
petus for the analysis presented here. 

1.2. Paper Outline 

This paper presents new secondary eclipse observations 
of WASP-12b's emission in a narrow band centered at 
2.315 /im, our detection and characterization of a cool 
star (which we call Bergfors-6) with high surface gravity 
near WASP-12, a correction of past eclipse measurements 
for the dilution caused by Bergfors-6, and our interpre- 
tation of WASP-12b's atmospheric emission. 



We describe our secondary eclipse observations and ini- 
tial data reduction in Sec. 2. As described in Sec. 3 we fit 
numerous model light curves to the data, select the sta- 
tistically optimal combination of parameters to use, and 
present the results of this eclipse. In Sec. 4 we describe 
our analysis of Bcrgfors-6's properties, and in Sec. 5 we 
use the results of this analysis to correct past transits 
and occultations of WASP-12b. In Sec. 6 we discuss our 
analysis of WASP-12b's corrected emission spectrum and 
provide updated constraints on the planet's bolometric 
luminosity. Finally, we conclude and suggest relevant 
possibilities for foUowup in Sec. 7. 

2. SUBARU/MOIRCS NARROWBAND TIME-SERIES 
PHOTOMETRY 

2.1. Summary of Observations 

We described recently the first tentative detection of 
emission from WASP-12b via spectroscopy at the 3 m 
NASA Infrared Telescope Facility (IRTF) (Crossfield 
et al. 2012). However, our precision was strongly lim- 
ited by chromatic and time-dependent slit losses result- 
ing from the use of a single, narrow (3") slit. We sub- 
sequently obtained time on the Multi-Object InfraRed 
Camera and Spectrograph (MOIRCS; Ichikawa ct al. 
2006; Suzuki et al. 2008) at Subaru Observatory to con- 
duct multi-object occultation spectroscopy of WASP- 
12b. A coolant leak at Subaru caused damage that pre- 
vented us from obtaining spectroscopy, so we operated 
the instrument in imaging mode using a custom nar- 
rowband filter. This filter (NB2315) is centered at ap- 
proximately 2.315 /im with a width at half maximum of 
27 nm^, and so is very well suited to probe the strong 
absorption feature predicted to lie at this wavelength by 
models used to infer a high C/0 ratio (Madhusudhan 
et al. 2011, their Figure 1). 

We observed one secondary eclipse of WASP-12b on 
14 Dec 2011 (UT). The start of observations was de- 
layed by instrument problems, but we managed to begin 
about half an hour before ingress and observed contin- 
uously thereafter. We observed at a position angle of 
330° and read out frames in correlated double sampling 
(CDS) mode with a constant integration time of 21 s per 
frame, using a readout speed of 8 and two dummy reads 
(to suppress a known, variable-bias effect; Katsuno et al. 
2003). These readout parameters result in substantial 
overhead penalties, and we averaged only one frame per 
61 s over our 6.5 hr of observations (which cover an air- 
mass range of 1.6—1.02—1.3). We recorded 388 frames in 
total. Conditions were nearly photometric, with stellar 
flux variations of 1-2% apparent. 

Following standard practices for high-precision pho- 
tometry (e.g., dc Mooij & Snellen 2009; Rogers et al. 
2009) we defocussed the telescope to spread the starlight 
over more pixels, thereby increasing observing efficiency 
and reducing the effect of residual flat fielding errors. 
The instrumental seeing improved throughout the night, 
and to avoid any substantially nonlinear detector re- 
sponse we added additional defocus to the telescope sev- 
eral times. Because the Subaru autoguider was inoper- 
ative we had to periodically apply manual offsets to the 
telescope tracking. The tracking was rather poor and de- 

^ A transmission profile of the NB2315 filter is available upon 
request from T.K. 
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Fig. 1. — Instrumental trends during our observations. From 
top to bottom: Chip 1 and 2 electronics box temperatures, rela- 
tive X and y motions, telescope focus encoder setting, WASP-12 
raw flux, WASP-12 relative flux, ratio of median sky background 
(pre-calibration) in Chips 1 and 2, and airmass. The dotted lines 
indicate the four points of contact corresponding to a circular orbit 
with our best-flt secondary eclipse center. The vertical dashed line 
corresponds to the onset of the anomalous trend apparent in the 
sky background and relative photometry: we exclude all data after 
this in our flnal analysis. 



spite our corrections we observed image drifts as large as 
1.2" (10 pixels); however, subsequent software develop- 
ment at Subaru has improved the tracking in the absence 
of the autoguider. The temperature of both detectors 
(as reported by the CHIPBDX FITS header keywords) in- 
creased from 76.2 K to a constant 77.0 K over the first 
1.5 — 2 hr. All these instrumental trends are shown in 
Figure 1, but wc ultimately find that they do not signif- 
icantly affect our photometry. 

2.2. Initial Data Reduction 

MOIRCS splits its field of view across two detectors, 
and we reduce the data from each detector indepen- 
dently. We calibrate the raw frames following the stan- 
dard MOIRCS reduction prescription, which proceeds as 
follows. MOIRCS returns the UT date and time at the 
beginning and end of each exposure in its FITS header. 
We convert these to B3T)tdb for our subsequent analysis 
(Eastman et al. 2010). We dark-subtract each frame and 
divide the result by the stack median of a set of dark- 
subtracted dome flats. Next, we correct our data for the 
intrinsic nonlinearity of infrared detector arrays (Vacca 



et al. 2004)^. Using a set of flat frames taken with typi- 
cal counts ranging from 2,000-22,000 ADU, we compute 
the median linearity correction coefficients (ai, . . . , a.4) 
of Vacca ct al. (2004)'s Eq. 20: 



(l + aix + a2X + aax + a4X 



(1) 



(where Cni is the ratio of an ideally linear signal to the 
measured signal) to be 

ai, 02, 03, 04 = (0.00347574, -0.00436064, -0.00111471, 0.00048908) 
for Chip 1, and 

ai, 02, as, 04 = (0.0063929, -0.0139614, 0.00548463, -0.00081818) 

for Chip 2. We define x as the measured ADU counts 
divided by 10^ to avoid very small coefficients. We then 
iteratively apply the correction algorithm in Vacca ct al.'s 
Eqs. 21-26, while further requiring that C„i is always 
> 1. Convergence typically occurs within 4-5 iterations. 
The nonlinearity correction is critical in our analysis: it 
changes our relative photometry by as much as 0.5% in 
some frames, and it slightly reduces our final, residual 
photometric RMS from 0.234% to 0.232%. 

At this point in the analysis substantial scattered 
and/or background light remains: we remove this by scal- 
ing and subtracting a median-combined set of median- 
normalized, dithered, dark-subtracted sky images. We 
follow these procedures independently for data from both 
channels; requiring an identical level of sky subtraction in 
both channels does not significantly change our results. 

We extract photometry using our own aperture pho- 
tometry package^, which uses bilinear interpolation to 
account for partial pixels while conserving flux. In each 
frame we extract subregions around each star, perform 
ID cross-correlations to measure relative stellar motions, 
and interpolate over hot pixels, stuck pixels, and any 
pixels more than 6a discrepant from their mean value. 
We then recenter the photometric apertures and perform 
standard aperture photometry. 

In imaging mode MOIRCS offers a roughly 4'x7' field 
of view split equally over two 2048^ IIAWAII-2 detec- 
tors, which allows several comparison stars to be fit 
into the WASP-12 field of view. Using stars more 
than about 1.8 mag fainter than WASP-12 decreases 
our final precision. Our large photometric apertures 
also require us to avoid choosing comparison stars with 
nearby companions. This leaves five comparison stars: 
2MASS stars 06302437+2937293 and 06303222+2937347 
(on Chip 1) and 06302377+2939118, 06301801+2939204, 
and 06302280+2938338 (WASP-12 is on Chip 2). Our fi- 
nal results are consistent (though of lower precision) if we 
use fewer comparison stars or use comparison stars falling 
only on a single detector. We examine the results from 
photometric apertures of various sizes and ultimately use 
target and inner and outer sky apertures with diameters 
of 39, 47, and 72 pixels (4.6", 5.5", and 8.4"). This choice 
minimizes the root mean square (RMS) of the residuals 
to our model fits; the final RMS is 0.232%. 

2.3. Instrumental Systematics 

^ A Python implementation of our MOIRCS nonlinearity cor- 
rection algorithm is available from the primary author's website. 
^ Available from the primary author's website. 
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We plot several variable instrumental parameters, 
along with the absolute and relative photometry of 
WASP-12, in Figure 1. One variable dominates in terms 
of its impact on our photometry: the curious trend in 
the ratio of the median sky background background mea- 
sured in the two detectors, which begins an anomalous 
excursion as WASP-12 crosses the meridian (only 10- 
15 min after egress) before later stabilizing at a new 
level. The relative photometry shows a qualitatively sim- 
ilar trend superimposed on a secondary eclipse (visible in 
the raw data). 

We also see this trend when dividing stellar photome- 
try from Chip 2 (excluding WASP-12) by photometry 
from Chip 1, and we even see it (at a lower ampli- 
tude) when comparing multiple reference stars on Chip 
2 against each other; it is thus a field-dependent effect. 
Because the trend begins just as WASP-12 crosses the 
meridian, we hypothesize that some loose component in 
the telescope or instrument settled in response to the 
change in the direction of the gravity vector. 

One possible culprit in this scenario is our narrow- 
band filter, whose spectral transmission profile depends 
on the angle of incidence of incoming light. To first 
order, increasing the angle of incidence translates the 
transmission profile to shorter wavelengths. The filter 
profile intersects a particularly strong telluric absorption 
(CH4) bandhead; from the vendor- supplied characteriza- 
tion data for our filter we estimate that a shift in the 
filter's angle of incidence of ~ 10° could induce a pho- 
tometric variation of the magnitude observed. However, 
no strong sky emission features are seen at these wave- 
lengths, so this scenario still has difficulty explaining the 
observed variation in the sky background. Regardless, 
subsequent MOIRCS multi-object spectroscopic data do 
not show this anomalous trend, a fact consistent with our 
hypothesis that the trend's presence is somehow related 
to the NB2315 filter. 

Whatever the cause of this anomalous trend, so long 
as we restrict our analysis to times before orbital phase 
0.5675 (the vertical dashed bar in Figure 1) our results 
change by less than 1.5(t no matter which comparison 
stars we choose. This choice leaves pre- and post-eclipse 
baselines which are rather short. We explored ways to 
use our entire data set by using the sky background trend 
as a decorrelation parameter (see Section 3 below), but 
such analyses resulted in larger fit residuals with substan- 
tially higher correlations on long timescales. We there- 
fore proceed by excluding the later data, while acknowl- 
edging the existence of this poorly-understood system- 
atic effect in MOIRCS narrowband imaging data. 

3. 2.315 fim NARROWBAND SECONDARY ECLIPSE 

We now present our analysis of the narrowband pho- 
tometry discussed in the preceding setion. In Section 3.1 
we describe the process of selecting an optimal model 
for our data and of fitting this model to the data. In 
Section 3.2 we describe the primary result of the fitting 
process: a 2.315 fim eclipse depth significantly discrepant 
from previous predictions. Then in Section 3.3 we show 
that the occultation we detect has a duration and time 
of center consistent with that expected for WASP-12b on 
a circular orbit. 

3.1. Fitting to the Data 



We fit our photometric time series with the following 
relation, representing a relative secondary eclipse light 
curve subjected to systematic effects: 

F,^ fo{l + dl,) (^ + J2c,v,}j (2) 

The symbols are: Fi, the relative flux measured at 
timestep i; /q, the true relative flux; ii, the flux in an 
occultation light curve scaled to equal zero out of eclipse 
and —1 inside eclipse; d, the normalized depth of sec- 
ondary eclipse; v^j, the J state vectors (i.e., image mo- 
tions, sky background, airmass, orbital phase, or low- 
order polynomials of these quantities) exhibiting a lin- 
early perturbative effect on the instrumental sensitivity; 
and Cj, the coefficients for each state vector. 

Experience shows that the choice of instrumental 
model is of crucial importance in extracting the most 
accurate system parameters from transit and occulta- 
tion observations (e.g., Campo et al. 2011). We therefore 
explore a large region of model parameter space by fit- 
ting our photometry using many different combinations 
of state vectors and a fixed secondary eclipse time and 
duration. We then use the Bayesian Information Crite- 
rion (BIC*) to choose which of these many models best 
represents our data. To do this we first assign uncertain- 
ties to each data point equal to the RMS of the residuals 
to an eclipse fit with no additional decorrelation vari- 
ables. We find that the BIC-minimizing model includes 
a transit light curve and a linear function of time, but 
no additional parameters. The model with the next-best 
BIC (3.2 units higher) also include a linear function of 
the X position of the stars on the detector; both this 
model and a model including no decorrelation parame- 
ters (ABIC ~ 10.7) return secondary eclipse depths con- 
sistent with that of our optimal model. 

We then again fit our preferred instrumental model to 
the data, but now allowing the secondary eclipse cen- 
ter, duration, and depth to vary while holding fixed the 
scaled semimajor axis (a/i?*) and the orbital inclination 
at the values listed in Table 1. We assess the uncertain- 
ties on the best-fit parameters using both the Markov 
Chain Monte Carlo and prayer bead (as described in 
Winn ct al. 2008, ; see also Jenkins et al. et al. 2002) 
approaches. The two sets of parameter distributions are 
quite consistent, which suggests correlated noise does not 
strongly affect our photometry. In both cases the re- 
sulting parameter distributions are unimodal, symmet- 
ric, approximately normal, and (excepting the standard 
correlation between occultation depth and baseline flux) 
uncorrelated. In the following we quote only the MCMC 
results, which provide substantially denser sampling of 
the posterior distributions than the prayer-bead results. 

3.2. Initial Narrowband Eclipse Depth 

We plot the data, our best-fit (linear baseline) model, 
and the residuals in Figure 2. The best-fit secondary 
eclipse depth is 0.41 % ± 0.05 %; note that the true 
planet/star flux ratio is ~ 10% greater than this, as 
we describe in Section 5 and Table 2. We find shal- 
lower or deeper best-fit eclipse depths (ranging from 

* Bayesian Information Criterion (BIC) = + klnN, where k 
is the number of free parameters and the number of data points. 
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TABLE 1 

WASP-12B: 2.315 /iM Narrowband Secondary Eclipse Parameters 



Parameter 



Units 



Value 



Reference 



P 
a/R, 

^offset 

758 
e costj 
esinoj 
Fp/F, (observed) 
Fp/Ff (corrected) 
Ts,2.3l5 



days 1.091423 
3.14 

BJDtdb 2455910.9090 ± 0.0013 
s 7±110s 
min 179.6 lb 4.5 min 

+0.00006 ± 0.00091 
0.011 ±0.013 
0.41 % ± 0.05 % 
0.45 % ± 0.06 % 
K 3640 K ± 230 K 



Hebb et al. (2009) 
Hebb et al. (2009) 
This work 
This work 
This work 
This work 
This work 
This work 
This work 
This work 



1.015 



3 1.010 

M— 

T3 
01 
N 



1.005 



1.000 



0.008 
I 0.004 
^ 0.000 





0.50 0.55 0.60 
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Fig. 2. — Top: relative 2.315 fim narrowband photometry of 
WASP- 12 (points, binned by a factor of ten for plotting purposes; 
errorbars are the standard deviation on the mean of each set of ten 
points) and our best-fit model (solid line). The Icr range of our 
model is also indicated by the dashed curves. Solid points are used 
in our analysis, while open points are excluded. Middle: Residuals 
to the fit. Bottom: differential sky background measured in the two 
MOIRCS detectors. The vertical dashed line indicates the onset of 
the photometric ramp apparent in sky and stellar photometry; wo 
exclude this data from our analysis, though the occultation depth 
is unchanged if we use all data and include the sky trend as an 
additional decorrelation parameter. The dotted lines indicate the 
four points of contact corresponding to a circular orbit with our 
best-fit eclipse center. 

roughly 0.41% to 0.50%) when using, respectively, larger 
or smaller photometric apertures, which indicates the 
sensitivity of this analysis to our particular choice of pa- 
rameters. For this reason we quote the measurement 
uncertainty of 0.05 % above, which is roughly twice that 
predicted from our prayer-bead analysis. 




13 10 30 

Number of Individual Points Binned 

ElG. 3. — RMS from binning the residuals shown in Figure 2 
over increasing numbers of data points (solid black line). This 
curve is exceeds the N~^^^ expectation from white noise (dashed 
line) by increasing amounts on successively longer timescales. The 
timescale of ingress or egress is indicated by the vertical dotted 
line. 

The residuals to the best fit have a relative RMS of 
0.232%. Photon (target -I- sky) noise considerations pre- 
dict a typical per-frame precision of 0.06%, so our per- 
formance is comparable to that obtained with other NIR 
secondary eclipse photometry (e.g., CroU et al. 2011). 
Figure 3 shows that our residuals bin down somewhat 
more slowly than N^^^^, indicating a moderate level of 
correlated residuals. 

The largest residuals occur at a "bump" in the light 
curve at orbital phase 0.47-0.48. One possible explana- 
tion is that this excursion is caused by telluric variations 
that are not entirely common mode across the MOIRCS 
field of view. An alternate explanation would be a short- 
term flare from Bergfors-6. We apply Difference Image 
Analysis (Bramich 2008)^ to two images generated by 
co-adding individual frames during and immediately af- 
ter the bump. The difference image shows no flux excess 
at the location of Bergfors-6 during this bump, indicating 
that the residual feature in the time series in not associ- 
ated with Bergfors-6. More exotic explanations, such as 
an attribution of this bump to accretion onto WASP- 12, 
should be treated with skepticism at present. 

" Our Python implementation of this algorithm is available at 
the primary author's website. 
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TABLE 2 

Dilution Factors and Corrected WASP-12b Transit and Occultation Depths 



Filter 



Reported Depth Aperture Fraction Dilution Fraction Corrected Depth Reference'' 



J 
H 

Ks 

Ks (MKO) 

NB2315 
IRAC CHI 
IRAC CH2 
IRAC CHS 
IRAC CH4 
IRAC CHI 
IRAC CH2 
IRAC CH2 

V/i"= 
IRAC CH2 
IRAC CHI 
IRAC CH2 
Johnson R 



0.00082 
0.00131 
0.00176 
0.00309 
0.00281 
0.0041 
0.00379 
0.00382 
0.00629 
0.00636 
0.0033 
0.0039 
0.0050 



± 0.00015 
± 0.00028 
± 0.00018 
± 0.00013 
± 0.00085 
± 0.0005 
± 0.00013 
± 0.00019 
± 0.00052 
± 0.00067 
± 0.0004 
± 0.0003 
± 0.0004 



0.01252 ± 0.00045 
0.0126 ± 0.0004 
0.0125 ±0.0003 
0.0112 ±0.0004 

0.01380 ± 0.00016 



1.000 
1.000 
1.000 
1.000 
1.000 
1.000 
0.902 
0.911 
0.855 
0.788 
0.850 
0.833 
0.833 



± 0.000 
± 0.000 
± 0.000 
± 0.000 
± 0.000 
± 0.000 
± 0.018 
± 0.016 
± 0.036 
± 0.068 
± 0.038 
± 0.049 
± 0.049 



0.230 ±0.100 
0.833 ± 0.049 
0.850 ± 0.038 
0.833 ± 0.049 
1.000 ±0.000 



0.0397 ± 
0.0606 ± 
0.0830 ± 
0.0981 ± 
0.0994 ± 
0.1002 ± 
0.1168 ± 
0.1204 ± 
0.1217 ± 
0.1307 ± 
0.1168 ± 
0.1204 ± 
0.1204 ± 



0.0017 
0.0029 
0.0039 
0.0047 
0.0047 
0.0045 
0.0055 
0.0060 
0.0060 
0.0063 
0.0055 
0.0060 
0.0060 



0.0186 ±0.0023 
0.1204 ± 0.0060 
0.1168 ±0.0055 
0.1204 ± 0.0060 
0.01571 ± 0.00096 



0.00085 
0.00139 
0.00191 
0.00339 
0.00309 
0.0045 
0.00419 
0.00424 
0.00694 
0.00701 
0.00363 
0.00429 
0.00550 



± 0.00016 
± 0.00030 
± 0.00020 
± 0.00014 
± 0.00093 
± 0.0006 
± 0.00014 
± 0.00021 
± 0.00057 
± 0.00074 
± 0.00044 
± 0.00033 
± 0.00044 



0.01257 ± 0.00045 
0.01386 ± 0.00044 
0.01374 ± 0.00033 
0.01232 ± 0.00044 
0.01402 ±0.00016 



LMIO 
Crll 
Crll 
Crll 
Z12 
This work 
Call 
Call 
Call 
Call 
Col2 
Col2 
Col2 (null)'' 

Cll<= 
Col2 (null)'' 
Col2 
Col2 
Mil 



B 
V 



0.0045 ±0.0015 
0.0090 ± 0.0030 
0.0281 ± 0.0014 



LMIO: Lopez-Morales et al. (2010), Crll: Croll ot al. (2011), Z12: Zhao et al. (2012), Call: Campo ct al. (2011), 
Col2: Cowan et al. (2012), Mil: Maciejewski et al. (2011), Cll: Chan et al. (2011) 

^ Results from the "null hypothesis" of Cowan ct al. (2012), which assumes zero ellipsoidal variation in their 4.5 /tm 
observations. 

^ These transit analyses average multiple photometric bands, so their correction factors may be less precise. Sec 
Sec. 5. 



We also attempted an alternative analysis in which we 
used all the data (including that affected by the anoma- 
lous background trend) and included the sky background 
trend as an additional decorrelation state vector. This 
analysis gives a marginally deeper secondary (0.53% ± 
0.05%), but has a higher residual RMS (0.261%), ex- 
hibits substantially higher levels of correlated noise when 
averaging on long time scales, and shows a significantly 
nonzero ecosw (inconsistent with previous analysis; cf. 
Campo et al. 2011; Croh et al. 2011; Cowan et al. 2012). 
These results suggest that our simpler model, which ex- 
cludes the latter part of our data set, gives the more 
reliable occultation measurement. 

3.3. Occultation Duration and Timing: No Surprises 

We find no evidence for significant deviations in sec- 
ondary eclipse duration or in the eclipse's time of cen- 
ter as compared to expectations from transit observa- 
tions and a circular orbit. We find a best-fit eclipse du- 
ration of 179.6 ± 4.5 min, and the eclipse occurs later 
than predicted by 7 ± 110 s (after accounting for the 
23 s light travel time across the system). Again, the pa- 
rameter distributions are unimodal and approximately 
normal. An analysis of previous Ks band observations 
reported a marginally longer secondary eclipse dura- 
tion (195 ± 7 min; Croll et al. 2011), while a weighted 
mean of Spitzer/IRAC occultations give a duration of 
177.7±2.1 min (Campo ct al. 2011, ; Cowan et al. do not 
report the durations of their transits and eclipses) and 
the z' occultation showed a duration of 169 min (with no 
uncertainty reported; Lopez-Morales ct al. 2010). Our 
narrowband measurement is consistent with these last 
two values (and with the duration expected from a cir- 
cular orbit) and is within 3a of the Ks band broadband 



results. Our data provide no evidence for an offset or 
longer-duration eclipse. 

Together, the secondary eclipse timing and duration 
tightly constrain the orbital eccentricity and longitude 
of periastron (Winn 2010; Scagcr 2011). We deter- 
mine ecosw and esinw to be +0.00006 ± 0.00091 and 
0.011 ± 0.013, respectively, which we interpret as being 
consistent with a circular orbit and with previous results 
based on eclipse and radial velocity observations (Campo 
et al. 2011; Croll et al. 2011; Husnoo et al. 2011). The 
time of eclipse also constrains the planetary velocity off- 
set expected at transit center, which can mimic wind- 
induced velocity offsets measured with high-resolution 
spectroscopy (Snellen et al. 2010; Fortney et al. 2010; 
Montalto et al. 2011; Miller-Ricci Kempton & Rauscher 
2011). Using Eq. 3 of Montaffo et al. (2011) we set a 3cr 
upper limit on any such orbit-induced velocity offset of 
+0.59 km s"\ 

4. Bergfors-6: AN OBJECT VERY CLOSE TO WASP-12 
4.1. Introducing Bergfors-6 

Before undertaking an analysis of WASP-12b's at- 
mospheric properties, we first pause to describe our 
improved characterization of a recently detected point 
source ~ 3 mag fainter than WASP-12 and only 1" away 
(Bergfors et al. 2011, 2012). This object requires us to 
revise upward past measurements of the planet's transits 
and occultations. 

During our Subaru observations one of us (I-T.) noticed 
a slight elongation in our (defocussed) images. This moti- 
vated us to refocus the system at the end of the night, and 
we recorded the image shown in Figure 4. It clearly shows 
a point source roughly 1" from WASP-12. A subsequent 
literature search revealed that this object was recently 
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TABLE 3 

WASP12/BERGFORS-6 ASTROMETRY 



32.90s 32.85s 32.80s 32.75E6h30m32.70s 



-1.0 -0.5 0.0 0.5 1.0 1.5 



RA(j20oo) A^[arcsec] 
Fig. 4. — Images used for astrometric and relative flux mea- 
surements: at left, seeing-limited image from Subaru/MOIRCS 
(2.315/im narrowband); at right, speckle image from IRTF/SpeX 
O^MKo)- Astrometric parameters derived from these images are 
listed in Table 3. Both images are displayed at the same orienta- 
tion and scale; they have different (logarithmic) color stretches in 
order to highlight the fainter companion. 

discovered using i and z photometry and assigned a pre- 
liminary spectral type of K4-M1 V (Bcrgfors et al. 2011, 
2012). We refer to this object as Bergfors-6, because it 
is the sixth object in Table 2 of Bergfors ct al. (2012). 
The existence of a bound companion at this projected 
separation (< 300 AU) would have potentially profound 
implications for the dynamical history of the system, and 
could provide a mechanism for Kozai-induced eccentric- 
ity and subsequent tidal heating to inflate the planet's 
radius to its present size (Fabrycky & Tremainc 2007; 
Nagasawa et al. 2008; Ibgui et al. 2011). 

Bergfors-6 has not been remarked upon in previous op- 
tical and infrared transit and occultation observations of 
the WASP-12 system (Hcbb et al. 2009; Lopez-Morales 
et al. 2010; Chan et al. 2011; Maciejewski et al. 2011; 
Campo et al. 2011; CroU et al. 2011; Cowan et al. 2012; 
Crossfield et al. 2012; Zhao et al. 2012). This is likely 
because with seeing-limited or Spitzer/IRAC resolution 
the two objects are at best only marginally resolved. The 
case is worse for most high-precision ground-based pho- 
tometry because of the common practice of substantially 
defocusing the telescope, which will clearly preclude de- 
tection of objects such as Bergfors-6. 

This star falls within the photometric apertures used 
in most previous analyses and dilutes the transit and 
secondary eclipse signals that have been measured (e.g. 
Daemgcn ct al. 2009). In this section we confirm the pre- 
vious detection of Bergfors-6 and more tightly constrain 
its spectral type. In the following section we then correct 
previous transit and occultation measurements for the 
photometric contamination of WASP-12 by Bergfors-6. 

4.2. Observations of Bergfors-6 
4.2.1. Subaru/MOIRCS NB231.5 Image 

As described above, we recorded the single well- 
focused MOIRCS image shown in Figure 4. We register 
the image's coordinate system using the 2MASS point 
source catalogue (Skrutskic et al. 2006) and confirm the 
MOIRCS plate scale to be 0.1 17 ±0.001" pix^^ (as listed 
in the instrument documentation). We conservatively 
adopt an uncertainty of 1° in the instrumental position 
angle. 

Bergfors-6 sits in the wings of the WASP-12 point 
spread function (PSF), so we perform a simultaneous 



Parameter 


Units 


Subaru/MOIRCS 


IRTF/SpeX 


Filter 




NB2315 


Ka/k o 


Flux ratio 




0.108 ± 0.007 


0.1048 ±0.0059 


Separation 


arcsec 


1.055" ± 0.026" 


1.078" ± 0.033" 


Position Angle 


deg 


250° ± 1° 


249.4° ± 1.1° 


Date 


UT 


2011 Dec 14 


2012 Feb 25 



fit to the two-dimensional PSFs of both stars. We use 
multiple elliptical Gaussian functions, holding the rota- 
tion and dispersion parameters fixed in each of the model 
PSFs and allowing only a single central location for each 
star. Thus for n Gaussian functions we have (5-|-5n) free 
parameters. We set the pixel uncertainties equal to the 
expectation from photon and read noise. 

We find that three elliptical Gaussians minimize the 
BIG, so we adopt this model and use Markov Chain 
Monte-Carlo (MGMG) techniques to explore the range 
of valid parameter space. We find the subsequent pa- 
rameter distributions to be unimodal and approximately 
Gaussian. To conservatively account for the uncertain- 
ties inherent in estimating accurate photometry and as- 
trometry from a single frame, we inflate the parameter 
uncertainties estimated from our MCMC by a factor of 
two. Our final determination of the flux ratio, separa- 
tion, and system position angle from the MOIRCS data 
are listed in Table 3. The astrometry is consistent with 
the initial discovery values (Bcrgfors et al. 2011, 2012). 
We confirm this 2.315 fim flux ratio by comparing aper- 
ture photometry of WASP-12 and (after subtraction of 
the best-flt WASP-12 PSF model) of Bergfors-6: this ap- 
proach gives a consistent result. The measurements pre- 
sented here are consistent with, but more precise than, 
estimates derived from several of our more poorly focused 
MOIRCS frames. 

4.2.2. IRTF/SpeX K Band Lucky Imaging 

The well-focused Subaru image described above mo- 
tivated us to acquire additional images of Bergfors-6. 
On 2012 Feb 25 (UT) we imaged the WASP-12 system 
with the IRTF/SpeX guide camera (Rayner et al. 2003), 
which uses a 512^ Aladdin 2 Insb array with a plate scale 
of 0.1185" pix~^. Observing through sometimes patchy 
clouds, we acquired 1,200 0.4 s Kmko (Tokunaga ct al. 
2002) frames (from airmass 1.2 - 1.5) and 900 0.21 s 
Jmko frames (from airmass 1.5—2.2). In all observations 
we held the SpeX instrument rotator at a position angle 
of 90°. The J band data were not sufficient to reliably 
detect Bergfors-6 (presumably because of the increased 
noise penalties resulting from the use of very short expo- 
sures and the smaller total integration time); hereafter 
we discuss only the K band data. 

We calibrate the SpeX images using a median stack 
of internal (thermal) fiat fields and subsequently per- 
form binlinear interpolation over a few noticeably bad 
pixels. We select the top 10% of all frames on the ba- 
sis of the peak pixel flux near the location of WASP-12, 
then use the "shift and add" algorithm to align and stack 
these frames (our S/N is too low for more advanced algo- 
rithms; Jcfferies & Christou 1993; Schoedcl et al. 2011). 
Changing the fraction of frames used in our analysis from 
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TABLE 4 

IRTF/SpeX Astrometric Calibrators^ 



WDS identifier 


Separation 


Position Angle 




[arcsec] 


[degrees] 


06295+3414 


4.25 + 0.14 


256.8 + 0.6 


06051+3016 


11.74 + 0.40 


177.3 + 0.6 


06508+2927 


6.60 + 0.22 


23.7 + 0.6 



^ All observations were made in the MKO K band on 



UT 2012 Fob 25. 

5% to 40% leaves our results unchanged within our esti- 
mated uncertainties. (However, our most selective anal- 
yses (which use only 1 — 2% of the data) show a hint of 
north-south elongation < 0.3". Though the resulting im- 
age is quite noisy, we note that the initial discovery image 
(Bergfors et al. 2011, 2012) also shows a similar elonga- 
tion. We recommend additional high-resolution imaging 
of Bergfors-6 to test this elongation.) The final image 
from our standard (10%) analysis is shown in Figure 4: 
in this image WASP-12 exhibits an axisymmctric PSF 
with a Strehl ratio of roughly 8% and a full width at half 
maximum of 0.33" (roughly a factor of 3 better than the 
seeing-limited resolution) . 

For astrometric reference we observed three known 
multiple systems taken from Version 2012-02-12 of 
the Washington Visual Double Star Catalog (WDS 
06295-1-3414, 06051+3016, and 06508+2927 Mason et al. 
2001) moderately near WASP-12, with comparable mag- 
nitudes to WASP-12, and with separations of 4-12". We 
took twenty 0.5-1 s frames of each system in the same 
region of the detector as our WASP-12 images, and in 
each frame we compute the centroids of both components 
using standard IRAF tasks. From these measurements 
and their dispersion we derive a SpeX guider plate scale 
of 0.116" ± 0.004" and an intrinsic field rotation (i.e., 
true position angle minus measured position angle) of 
—0.5° ± 0.6°. We adopt these values in our subsequent 
analysis and list our astrometry of the WDS stars in Ta- 
ble 4. 

We determine the fiux ratio of the two stars using aper- 
ture photometry. Bergfors-6 is located in the PSF wings 
of WASP-12, so we must account for this contamination. 
Because our PSF is quite symmetric (though distinctly 
non-Gaussian) we compute an average radial profile for 
WASP-12 (after masking out the 90° wedge of sky di- 
rected toward Bergfors-6). We reinterpolate this one- 
dimensional profile into a two-dimensional model PSF. 
We estimate the uncertainty of the profile by taking the 
standard deviation on the mean in each annular bin, and 
propagate these uncertainties along with the combined 
photon and read noise. We then subtract the WASP- 
12 model PSF from the image and compute partial-pixel 
aperture photometry at the locations of Bergfors-6 and 
WASP-12. Residuals are still apparent near the center 
of WASP-12, so we restrict our analysis to smaller aper- 
tures: an inner aperture radius of 2.5 pix provides the 
highest S/N (and least evidence for contamination) for 
Bergfors-6, so we use this aperture for both systems. Our 
final estimate of the Kmko flux ratio is listed in Table 3, 
and it is consistent with our narrowband MOIRCS mea- 
surement. 

We measure the relative astrometry of WASP-12 and 



Bergfors-6 by computing the centroid of WASP-12 in the 
speckle image, and of Bergfors-6 in the profile-subtracted 
image. We estimate the uncertainties in these measure- 
ments by bootstrap resampling (Press 2002), in which we 
repeat our analysis many times using synthetic data sets, 
constructed by sampling (with replacement) our original 
set of 1,200 images. We list the separation and position 
angle derived from the IRTF speckle data in Table 3. 

4.2.3. Spitzer/IRAC Imaging 

We also examined Spitzer/IRAC subarray data 
(3.6 and 4.5 fj,m, from Cowan et al. 2012) to search for 
evidence of Bergfors-6. We performed a weighted least 
squares fit to each median stack of 64 subarray frames 
(using the pixel uncertainties provided by the IRAC cal- 
ibration pipeline. Version 18.18.0) by linearly interpo- 
lating the appropriate 5x oversampled point response 
functions^" (PRF) to account for subpixel motions. 

We do see evidence for an additional point source in 
the IRAC data, located approximately 1-2 pixels west- 
southwest of WASP-12. However, we are unable to mea- 
sure precise astrometry or relative photometry with these 
data for several reasons. First, the IRAC plate scale 
(1.09" pix-i) is comparable to the WASP-12/Bergfors-6 
separation; second, the IRAC PSF is undersampled at 
these wavelengths. Consequently, we see clear evidence 
for oversubtraction in the PRF fitting at the location of 
Bergfors-6, so we cannot reliably determine the system 
flux ratio (de-weighting the pixels closest to Bergfors-6, 
but offset from the WASP-12 core, does not change this 
result). From these measurements we estimate a fiux ra- 
tio of > 7% in the two IRAC channels, consistent with 
our ultimate interpretation of Bergfors-6 as a cool stellar 
object. 

4.2.4. Keck/NIRSPEC Spectroscopy 

We searched online data archives for additional 
evidence of Bergfors-6, and found a set of high- 
resolution K band spectra taken with Keck/NIRSPEC, 
a high-resolution, cryogenic, echelle, NIR spectrograph 
(McLean ct al. 1998), on UT 2010 Apr 22 (Keck Pro- 
gram ID C269NS, P.I. G. Blake). This data set con- 
sists of 16 four-minute integrations of WASP-12 taken 
using the 0.432" x 24" slit. The WASP-12 observations 
were taken at a position angle of ^73° (roughly aligned 
with WASP-12 and Bergfors-6), and the seeing was suffi- 
ciently good to distinctly resolve the two components in 
the spectra. 

We extract our spectra using our own set of Python 
tools to trace the spectra in the dark-corrected and fiat- 
fielded NIRSPEC frames. In each echelle order of each 
frame, we compute a high S/N mean spectral profile by 
collapsing the trace along the dispersion direction and 
fit two Gaussian functions to this profile: this provides 
an estimate of the projected separation of WASP-12 and 
Bergfors-6 in each frame. We then fit two Gaussian func- 
tions to each resolution element while holding constant 
the positions of the two sources: the amplitude of each 
Gaussian represents the flux in that wavelength element. 
We then compute weighted means from the individual 

1° Available at http://lrsa.ipac.caltech.edu/data/SPITZER/ 
docs/irac/calibrationf iles/psf prf / 
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extracted spectra and estimate uncertainties by measur- 
ing the variations in each pixel, after excluding points 
deviating by > 3a. 

Using a high-resolution simulated telluric spectrum 
(generated using ATRAN; Lord 1992) we identify known 
telluric lines and compute a best-fit dispersion function 
in each echelle order. Estimating a line centroid preci- 
sion of 0.5 pix, we find a cubic or quartic polynomial 
minimizes the BIC of these fits. The residuals to our dis- 
persion solutions have RMS values ^0.1 Aand maximum 
excursions of < 0.2 A. 

In the raw NIRSPEC frames the spatial axis of the slit 
is not aligned with the NIRSPEC detector columns, so 
spectra taken at the A and B nod positions are offset 
from each other. We spline-interpolate the spectrum in 
each echelle order and cross-correlate it at sub-pixel in- 
crements with a high signal to noise (S/N) template spec- 
trum (Doming et al. 2005). We construct our template 
by taking the temporal average, after removing outliers, 
of all our spectra. A parabolic fit to the peak of each 
spectrum's cross correlation provides the optimal offset 
value, and we then spline-interpolate all the spectra to a 
single, common reference frame. We then combine the re- 
sulting set of aligned spectra (excluding outlying points) 
and thereby provide a set of simultaneous high-resolution 
spectra of both WASP-12 and Bergfors-6. The spectra 
of WASP-12 and Bergfors-6 have median S/N values of 
203 and 32 per pixel, respectively. 

Because the two spectra are obtained simultaneously 
and the objects are separated by only 1" we expect the 
telluric signature in both spectra to be indistinguishable. 
We therefore divide the spectrum of Bergfors-6 by that 
of WASP-12 to remove the effect of telluric absorption. 
Possible misalignment of the spectrograph slit prevent 
these data from usefully constraining the absolute flux 
ratio of these two objects, but the data constrain tightly 
Bergfors-6's spectral type from the relative strengths of 
individual spectral features. The NIRSPEC spectra do 
not cover the wavelengths of standard gravity indicators 
such as Na or Ca lines, so our subsequent analysis focuses 
on the two most prominent gravity-sensitive features cov- 
ered by these data: the ^^CO (2, 0) and (4, 2) bandheads 
located at 2.294 /j,m and 2.353 /xm (Kleinmann & Hall 
1986). As we now describe, we find that Bergfors-6 is a 
hot M dwarf. 

4.3. The Spectral Type and Nature of Bergfors-6 
4.3.1. Photometric Constraints 

With our four (i, z, KmkOi 1^2315) relative photomet- 
ric measurements and the NIRSPEC spectrum we de- 
termine the spectral type of Bergfors-6, as described be- 
low. First, we apply the relationship between spectral 
type and absolute magnitude of Kraus & Hillenbrand 
(2007) using our photometry. This relationship is for 
main sequence stars, and WASP-12 is 25% larger than 
a zero-age main sequence star of the same mass (Hebb 
et al. 2009; Torres et al. 2010). Accounting for this, 
and assuming Te/f — 6300 ± 100 K, gives a distance 
modulus for WASP-12 of yLtwi2 = 7.7 ± 0.2 mag (signif- 
icantly nearer than the previous estimate; Chan et al. 
2011). For Bergfors-6, our estimate of the i — K color 
(after applying the color transformations of Carpenter 
2001) implies a main-sequence spectral type of MO-Ml 



(Te// = 3700 ± 100 K). Assuming systematic uncertain- 
ties of 0.2 mag gives = 7.1 ± 0.2 mag, rather closer 
than WASP-12 if Bergfors-6 is still on the main sequence. 
That Bergfors-6 is closer to Earth than is WASP-12 is the 
opposite of the trend noted by Dacmgen ct al. (2009), 
who uniformly estimated that their faint companions to 
planet host stars were more distant than the brighter 
component. 

Next, we fit only the relative (Bergfors-6/WASP-12) 
stellar photometry using low-resolution stellar atmo- 
sphere models (Castelli & Kurucz 2004). We interpolate 
to the effective temperature, [Fe/H], and surface gravity 
of WASP-12 and hold these values fixed in the model- 
ing. For Bergfors-6 we assume a metallicity equal to 
that of WASP-12 and allow three free parameters: sur- 
face gravity, effective temperature, and a geometric fac- 
tor (f =-MB6_dYii2 \ relating the relative sizes and helio- 

centric distances of the two stars. A standard Pythonic 
minimizer and an MCMC analysis using the emcee affine- 
invariant sampler (Forcman-Mackcy ct al. 2012) provide 
the desired physical parameters and their uncertainties. 
The derived parameters for Bergfors-6 are 3840 ± 70 K 
(which agrees well with our previous estimate of this ob- 
ject's effective temperature) and / =0.452 ± 0.015; the 
covariance between these parameters is -1.73 K. 

We also considered that Bergfors-6 might be an ex- 
tragalactic, rather than a stellar, contaminant (Luhman 
& Mamajek 2010). However, a comparison of its pho- 
tometric spectral energy distribution with low-resolution 
galactic spectral templates (Assef et al. 2010) suggests 
that this explanation is unlikely. 

4.3.2. Spectroscopic Constraints 

Spectroscopy is particularly well-suited to constrain 
surface gravity. We now use both the NIRSPEC spec- 
trum and the relative photometry described above to 
constrain Bergfors-6's parameters, again using the emcee 
MCMC sampler (Foreman-Mackey et al. 2012). For this 
analysis we use the BT-Settl library^ ^ computed using 
the PHOENIX atmosphere code(Allard et al. 2010). This 
library provides high-resolution model spectra across a 
wide range of parameter space. We use the so-called 
"hot" models using abundances from Asplund et al. 
(2009) with no alpha enhancement. As noted previously, 
we use only the two NIRSPEC echelle orders that cover 
CO bandhead features - the wavelengths from 2.273- 
2.308 ^m and 2.344-2.380 ^m. 

For a given set of input parameters, our spectral mod- 
eling algorithm begins by logarithmically interpolating 
between BT-Settl models at the nearest values of Te//, 
log 5, and [M/H]. The model then applies (a) a Doppler 
shift, (b) a quadratic continuum normalization (because 
the absolute slope and curvature of the spectrum is un- 
known owing to possible slit misalignments), and (c) a 
convolution with a Gaussian kernel of specified width. 
Finally, we bin (not interpolate) the model spectrum onto 
our NIRSPEC pixel grid and compute relative broadband 
photometry as described in the preceding section - ex- 
cept that here we propagate the uncertainties in WASP- 
12's parameters into the modeling by performing a ran- 
dom draw from normal distributions in Tg//, logg, and 

Available online at http://phoenix.ens-lyon.fr/ 
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Fig. 5.— Keck/NIRSPEC spectrum of Bergfors-6 (black 
curve, top and middle) and relative photometry of WASP-12 and 
Bergfors-6 (filled points, bottom), and our best-fit Phoenix/BT- 
Settl model (red curve). As discussed in Section 4.3, the ensemble 
of measurements suggests Bergfors-6 is a hot M dwarf located 50% 
closer to Earth than WASP-12. Top and Middle: The blue curves 
show our estimated spectroscopic measurement uncertainties; the 
vertical scale for these curves is indicated at right. Bottom: The 
open points at bottom are the inferred photometric dilutions of 
transits or occultations measured in various bandpasses (tabulated 
in Table 2); we indicate filters used in this analysis with solid lines, 
while other filters are denoted with dashed lines. The error bars of 
the open points represent the 68.3% confidence intervals on these 
dilution estimates, taking into account the uncertainties in our fit. 

[M/H] at each step in the MCMC analysis that foUows. 

In our analysis the spectral type of Bergfors-6 is con- 
strained almost entirely by the ^ 2000 spectroscopic data 
points, while the geometric ratio 7^^^ is constrained 
only by the broadband photometry. The results of this 
analysis for Bergfors-6 are an effective temperature of 
3660^60 K and a \ogg (cgs) of 5.13iQ;22; we show our 
best-fit model spectrum in Figure 5. Our derived pa- 
rameters are fully consistent with an MO dwarf on the 
main sequence, which would imply a radius of 0.5— O.6i?0 
(Torres et al. 2010). The geometric ratio from our analy- 
sis is 0.520^Q 037, which implies a radius 50% larger than 
expected for a main-sequence dwarf lying at the same 
distance at WASP-12. Our spectral analysis therefore 
suggests that Bergfors-6 lies approximately 50% closer 
to Earth than does WASP-12 and that it represents a 
chance foreground alignment. 

4.4. Radial Velocities 

To estimate the radial velocities of the two stellar com- 
ponents in our telluric-corrected Bergfors-6/WASP-12 
ratio spectrum we cross-correlate with BT-Settl mod- 
els (Allard et al. 2010). We cross-correlate each echelle 
order of the ratio spectrum with models with effective 
temperatures of 6200 K and 3800 K. We also cross- 
correlate the spectra (before telluric correction) with the 
high-resolution atmospheric transmission profile of Hin- 
kle et al. (2003) to establish our observational reference 
frame. 

The radial velocity of Bergfors-6 is constrained mainly 
by the strong CO bands lying redward of 2.29 /im, while 
WASP-12's is tightly constrained only by the broad 



(FWHM = 8.lA) Brackett 7 hne. After correcting for the 
Earth's velocity along the line of sight (using the Python 
routine astrolib.baryvel) we estimate radial velocities 
for WASP-12 and Bergfors-6 of 16.5 ± 2.6 km s~^ and 
19.7 ± 1.3 km s~^, respectively. These values are consis- 
tent with the radial velocity of 19.1 km derived from 
WASP-12b's initial radial velocity measurements (Hebb 
et al. 2009; Campo et al. 2011). This common veloc- 
ity is consistent with a scenario in which either WASP- 
12 and Bergfors-6 are gravitationally bound and share a 
common three-dimensional space motion, or in which the 
consistency of the two stars' radial velocities is merely a 
coincidence. 

We also find no evidence for multiple spectral line pro- 
files, which would indicate Bergfors-6 is an unresolved 
binary. The cross-correlation profiles of our data and 
spectral template have full-widths at half-maximum of 
approximately 14 km s~^. All three echelle orders con- 
taining strong CO features (centered on 2.29, 2.36, and 
2.44 //m) show unimodal cross-correlation peaks and no 
evidence of the secondary peaks that would suggest an 
additional cool companion. Thus Bergfors-6 shows no 
spectroscopic evidence of binarity. 

4.5. Interpretation of Bergfors-6 

We conclude that Bergfors-6 is a cool star showing ab- 
sorption features consistent with a high surface gravity. 
Our spectroscopy implies that Bergfors-6 is a main se- 
quence star, and our relative photometry suggests that 
Bergfors-6 hes 50% closer to Earth than does WASP- 
12. In this scenario, the hint of elongation alluded to in 
Sec. 4.2.2 is spurious and the consistent radial velocities 
of WASP-12 and Bergfors-6 is coincidental. However, if 
Bergfors-6 were a binary M dwarf system observed near 
conjunction then the binary and WASP-12 could lie at 
the same distance from Earth; if the three components 
were bound this scenario would offer a natural explana- 
tion for the consistent systemic velocities. Because of 
the intriguing possibilities inherent in a gravitationally 
bound arrangement, we briefly discuss the implications 
of such a scenario below. 

At the distance of WASP-12, the projected separation 
of Bergfors-6 is roughly 400 ± 100 AU. Such an object 
would have an orbital period of several thousand years, 
and as such would be marginally compatible with the 
inferred upper limit of binary separation {'^ 300 AU) 
needed to substantially influence planetary migration 
and dynamics (Desidera & Barbicri 2007). Assuming 
Bergfors-6 is near apastron implies a Kozai oscillation 
timescale (Fabrycky & Tremaine 2007) of (very roughly) 

20[(1 - eB6)/(l + eee)]^/^ Gyr, where eee is Bergfors-6's 
orbital eccentricity. For substantial Kozai interactions to 
have taken place during WASP-12's lifetime (Hebb et al. 
2009; Chan et al. 2011) we thus require eee > 0.7. Long- 
period binaries with such high eccentricities are rare, but 
they do exist (Duquennoy & Mayor 1991). 

Fabrycky & Tremaine (2007) predict that hot Jupiter 
systems with an additional, widely-separated stellar com- 
ponent will preferentially exhibit misalignment between 
their stellar spin and planetary orbital axes. Recent ob- 
servations of WASP-12 have determined that the sky- 
projected angle between the spin and orbital axes is 
59^l^deg, strongly suggesting a misaligned system (Al- 
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brecht ct al. 2012). This may be circumstantial evidence 
that WASP-12 has migrated via Kozai interactions and 
that a long-period bound companion is required (e.g., 
Ibgui ct al. 2011). 

Ultimately, high-resolution imaging can most quickly 
determine whether WASP-12 and Bergfors-6 truly ex- 
hibit common proper motion and are gravitationally 
bound, and whether Bergfors-6 is a single or multiple sys- 
tem. Based on WASP-12's proper motion 8 mas yr^^; 
Zacharias et al. 2004), speckle or seeing-limited astrom- 
etry of the type presented here will not be sufficient for 
this purpose. However, a two-year baseline of large- 
aperture adaptive optics imaging (Yelda ct al. 2010) 
could suffice to confirm or rule out common proper mo- 
tion. 

5. REVISING PAST TRANSIT AND SECONDARY ECLIPSES 

Because Bergfors-6 was not noted in previous transit 
and secondary eclipse observations of WASP-12b, these 
flux diminutions were diluted by this faint star's con- 
stant baseline flux. This effect is largest in the infrared; 
the results of optical observations change by only a few 
percent, less than their typical uncertainties. Although 
a full re-evaluation of WASP-12's system parameters is 
beyond the scope of this work, we correct the depth mea- 
surements for the contamination effect and present re- 
vised transit and occultation depths below. We propa- 
gate the uncertainties in Bergfors-6's effective tempera- 
ture into our estimates of the photometric dilution caused 
by Bergfors-6, which we list in Table 2. In some sec- 
ondary eclipses, and in all transits, the corrections we 
apply change the previously reported depths by > la. 

Essentially all the light from Bergfors-6 lies within 
the apertures of ground-based observations (Hcbb ct al. 
2009; Lopez-Morales ct al. 2010; CroU ct al. 2011; Chan 
et al. 2011; Maciejewski et al. 2011; Zhao et al. 2012), but 
Spitzer/IRAC analyses use narrower apertures (Campo 
et al. 2011; Cowan et al. 2012) and so only a portion 
of Bergfors-6's starlight contaminates these secondary 
eclipse measurement. To estimate the IRAC contami- 
nation fraction we generate 10 x super-sampled PSFs for 
all four IRAC channels^^ , using a 6300 K blackbody spec- 
trum simulated at the center of the instrument field of 
view. We then compute aperture photometry using the 
reported photometric aperture diameters (Campo et al. 
2011; Cowan et al. 2012) at a position offset by 1.05" from 
the PSF center to estimate how much of Bergfors-6's flux 
fell into the WASP-12 aperture in these analyses. We ig- 
nore possible time-variable illumination caused by the 
intrapixel effect (Charbonneau ct al. 2005). Note that 
the phase curve observations of Cowan et al. (2012) are 
also diluted by Bergfors-6, and must be revised upward 
by the same factors as indicated in Table 2. This in 
turn increases the ellipsoidal variation inferred from these 
measurements, placing the IRAC 4.5 fim results in even 
stronger conflict with those from WFC3 (Swain ct al. 
2012). 

The correction of optical transit measurements is com- 
plicated by the common, but deplorable, practice of re- 
porting a single transit depth when using observations 
taken in different bandpasses (Hebb et al. 2009; Chan 

Using Tiny Tim; available at http : //ssc . spitzer . caltech. 

edu/ 



et al. 2011). Atmospheric characterization via transit 
observations depends on the fundamentally wavelength- 
dependent planetary radius during transit, and so we 
recommend that future analyses of multiple photomet- 
ric data sets report the transit depths measured in each 
bandpass in addition to a single, achromatic value. In 
addition, such analyses are not always clear about the 
relative weighting of data points from separate observa- 
tions. Nonetheless we attempt to estimate the relative 
weightings and derive appropriate correction factors for 
prior multi-band optical transit observations of WASP- 
12. 

The analysis of Maciejewski et al. (2011) used only 
a single bandpass (Johnson R) and so their Johnson R 
transit depth is the most reliable optical transit measure- 
ment in Table 2. The analysis of Chan ct al. (2011) uses 
two transit data sets: 671 V band and 470 i' band obser- 
vations with residual RMS values of 2.0 and 1.2 mmag, 
respectively. This work uses an additional multiplicative 
term to increase the per-point data uncertainties (1.48 
and 1.57, respectively). Assuming that the statistics of 
the final reported transit depth behaves similarly to a 
weighted mean, we estimate that the two data sets con- 
strain the transit parameters with roughly equal weight. 

We use these weights to determine a weighted aver- 
age of the correction factors for each bandpass. Such a 
weighted average is only a rough approximation to the 
true correction factor, so the correction factors for these 
analyses are somewhat less certain. We thank the referee 
for pointing out that the situation is even more muddled 
for the discovery paper (Hebb et al. 2009), which uses 
three transit data sets: 227 B band, 614 z' band, and 
6393 SuperWASP (roughly V band) observations. Be- 
cause Hcbb ct al. (2009) do not report the residual RMS 
scatter of the SuperWASP photometry, we do not at- 
tempt to determine the relative weighting of these sev- 
eral transit observations and we do not report a mean 
dilution-corrected transit depth for these data. 

The dilution correction factors, and dilution-corrected 
transit and occultation depths, are listed in Table 2. The 
optical transit depths in particular are increased by a 
few percent, and the planetary radius increases by only 
half this factor. Owing to current uncertainties in stellar 
properties, such an effect is smaller than current uncer- 
tainties on WASP-12b's physical radius. 

6. EMISSION SPECTRUM AND ATMOSPHERIC 
PROPERTIES 

We now return out attention to the nature of WASP- 
12b's atmosphere as constrained by its spectral energy 
distribution (SED). In our narrow bandpass we find a 
brightness temperature of 3640 K ± 230 K by using the 
model stellar spectrum described in the following sec- 
tion, modehng WASP-12b's emission in our bandpass as 
a blackbody, and propagating the la uncertainties in our 
occultation measurement. This brightness temperature 
is rather higher than the planet's equilibrium tempera- 
ture of 2990 ± 110 K (computed using the Bond albedo 
and recirculation efficiencies from Cowan ct al. 2012) and 
is higher than inferred in any broad photometric band- 
pass (Lopez-Morales et al. 2010; CroU ct al. 2011; Campo 
ct al. 2011; Cowan et al. 2012). The CFHT/WIRCam 
Ks filter cuts off just where the NB2315 filter cuts in, 
so our measurement is not in conflict with this previous 
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secondary eclipse observation (CroU ct al. 2011). 

In the following, we use the weighted averages of the 
two sets of IRAC (3.6 /xm and 4.5 /xm) measurements 
(Campo et al. 2011; Cowan et al. 2012, ; see also Table 2). 
Two possible 4.5 /xm secondary eclipse depths were re- 
ported by Cowan ct al. (2012), and we use their "null 
hypothesis" consistent with no ellipsoidal variations (as 
implied by the constraints placed on WASP-12b's shape 
by Swain et al. 2012). 

We first introduce the improved, model-independent 
constraints these observations place on WASP-12b's 
bolometric luminosity in Sec. 6.1, and we then discuss 
our efforts to generate a coherent model of the planet's 
emission spectrum in Sec. 6.2. 

6.1. Bolometric Luminosity 

WASP-12b has one of the best-determined bolomet- 
ric luminosities of any extrasolar planet (Cowan & Agol 
2011). In an earlier work (Crossfield et al. 2012) we dis- 
cussed the current constraints on the bolometric lumi- 
nosity of WASP-12b. Here we update this analysis in 
light of our dilution-corrected occultation measurements 
and recent Spitzer/IRAC and Hubble/WFCS observa- 
tions (Cowan et al. 2012; Swain et al. 2012). 

Current measurements now constrain the planet's day- 
side to have Lboi = (3.6 — 5.0) x lO'^^erg s~^, where the 
lower limit assumes the case of zero emission between 
the observed bandpasses. Here we have followed the 
same approach as our previous calculation, but we now 
use the WFC3 spectrum (Swain et al. 2012) from 1.1- 
1.485 ^m (replacing the J band eclipse, but retaining the 
H band eclipse, of CroU et al. 2011). In this analysis our 
narrowband measurement, and the controversy over the 
IRAC 4.5 /im measurement, affect the luminosity by only 
0.1 X lO^^erg s-^ 

On the basis of thermal occultation and phase curve 
measurements, WASP-12b's bolometric albedo has been 
inferred to be Ab = 0.25±0.1 (Cowan et al. 2012), which 
imphcs that the planet absorbs (3.8 ±0.8) x lO'^^erg s~^ 
from its host star. Following the approach of Crossfield 
et al. (2012), we find that this value constrains the night- 
side luminosity to be < 1.6 x 10^°erg s""'^, consistent with 
the nightside luminosity inferred by Cowan ct al. (2012) 
of 0.06ig j2 X lO^Oerg 8"^ This last point further as- 
sumes that the night side (like the day side; see below) 
emits approximately like a blackbody. 

We are thus nearing a bolometric luminosity suffi- 
ciently well-constrained that we can test evolutionary 
models of this planet. The current uncertainty in Lboi is 
dominated by measurements at the shortest wavelengths, 
suggesting that occultation observations at wavelengths 
< 1 /xm may be the best next step toward an even more 
tightly constrained bolometric luminosity. 

6.2. Atmospheric Models 

The first systematic effort to retrieve WASP-12b's at- 
mospheric parameters (using infrared broadband sec- 
ondary eclipse photometry) inferred a high C/0 ratio 
(> 1) and ruled out any strong temperature inversion 
at the pressures probed (0.01-2 bar; Madhusudhan ct al. 
2011). That study published several representative mod- 
els, all of which had 10 with seven measurements 
and ~ 10 free parameters (BIC~ 32). These models also 



predicted a strong absorption feature (depth < 0.2%) at 
2.315 /im: our narrowband secondary eclipse measure- 
ment rules out this absorption feature at > 3cr. This 
narrowband result, our characterization of Bergfors-6 
and correction for its diluting effect, and the alterna- 
tive Spitzer/IRAC 4.5 //m eclipse depth of Cowan et al. 
(2012) all indicate the need for a new analysis of WASP- 
12b's atmospheric properties. 

To better understand the nature of the planet's dayside 
emission we constructed a variety of atmosphere mod- 
els for WASP-12b, following Barman ct al. (2001, 2005). 
We compared model SEDs to the corrected occultation 
depths from Table 2. For this exercise we adopt the 
Cowan et al. (2012) 4.5 /xm value without ellipsoidal vari- 
ations (their "null hypothesis"). As discussed by Cowan 
et al., when ellipsoidal variations are allowed, the inferred 
planet elongation in this band is substantially larger than 
expected. Indeed, recent NIR HST/WFC3 observations 
exclude the large ellipsoidal variations inferred from the 
Spitzer/IRAC analysis (Swain et al. 2012). Further- 
more, the measured primary transit radius at 4.5 /im 
differs dramatically from radii at other wavelengths as 
well as predicted transit spectra. Given that two very 
unexpected physical properties are required to explain 
the ellipsoidal variations, the 4.5 /xm secondary eclipse 
value without ellipsoidal variations is likely more reliable 
(Cowan 2012, private communication). Finally, our ul- 
timate result - that the high C/0 result is not justified 
by current photometric data - does not change even if 
we use the average of the Campo ct al. and Cowan et al. 
results. 

By far the model that best matches the broad band 
occultation photometry is a 3000 K black body (upper 
panel of Fig. 6, and Fig. 7), with = 15 with only 
two free parameters (BIC= 19). Taking the mean of the 
Campo et al. and Cowan et al. results gives a black body 
model with higher = 25 and BIC= 29, but this BIC 
value is still smaller than that for the high C/0 models. 
Discussed above. This lower BIC value indicates that 
current observations do not justify the use of these more 
complicated models. Our best-fit black body temper- 
ature matches the expected equilibrium temperature of 
the planet if redistribution of heat to the night side is ex- 
tremely inefficient. The close match to a blackbody also 
indicates that the planet's photosphere is nearly isother- 
mal at the depths probed by these wavelengths. 

Our default irradiated atmosphere model assumes zero 
redistribution of heat to the nightside (consistent with 
the black body analysis) and metal abundances matching 
those of the host star ([Fe/H] = 0.3). This model pre- 
dicts a temperature inversion above an isothermal region 
(Fig. 7); it also has = 15; since it has substantially 
more free parameters than a black body, the BIC for 
this model fit is substantially worse than for the simpler 
model in the previous paragraph. In such a model the 
inversion extends partially across the IR photosphere re- 
sulting in modestly inverted CO and H2O bands (lower 
panel of Fig. 6). However, the model band-integrated 
fluxes are similar to the black body and match the broad 
band observations equally well (Fig. 6) . Naturally, then, 
the blackbody model provides a substantially lower BIC 
than the radiative-transfer models (both our model and 
the models presented by Madhusudhan et al. 2011). 
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2 4 6 S 10 

Fig. 6. — WASP-12b emisaoii ' spi'ct^jim^; see Sec. 6 for a full 
discussion. Solid points are the dilution-corrected photometric sec- 
ondary eclipse depths listed in Table 2; at 4.5 fim we plot the results 
of both Campo ct al. (2011) and Cowan ot al. (2012) The lower 
panel is the black body comparison (open symbols are the model, 
band-integrated, points). The upper panel includes several spec- 
tra: the dark black spectrum (with open black symbols) is our solar 
abundance, 2tt redistribution model spectrum (smoothed for plot- 
ting purposes by convolving with a Gaussian of FWHM = lOOA). 
The grey spectrum reproduces the black model but is plotted at 
100 times higher resolution to show the many narrow emission lines 
predicted by this model. The inset shows the narrow band region 
and a model with CO lines scaled up (red curve). Within the inset 
plot, the "standard" high-res model is also shown in the same grey 
color as in the main figure. Note the scale: these lines are narrow 
but very bright. If such strong, narrow emission lines are present, 
they should be easily discerned by future observations. 
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Fig. 7. — Temperature-pressure profile (top) and molecular abun- 
dances (bottom) used to construct the model emission spectrum of 
WASP-12b shown in Fig. 6. The middle panel shows the normal- 
ized contribution functions for the indicated filters. 

Although our narrowband measurement is fit by the 



simple blackbody model (within 2a) , we additionally in- 
vestigated what mechanisms might cause such a narrow 
band to exhibit a substantially higher flux. In an attempt 
to reproduce the narrow band flux, we made a number 
of ad-hoc modifications to the temperature-pressure (T- 
P) profile of the default atmosphere model. The base of 
the inversion was moved to higher and lower pressures, 
the steepness of the inversion was altered and several 
Gaussian temperature perturbations were added. None 
of these models could match the narrow-band flux with- 
out negatively impacting the comparison at other wave- 
lengths. A dozen different carbon and oxygen abun- 
dances (ranging from 7.7 — 9.4 on the standard base- 
10 logarithmic scale, where H2 abundance is exactly 12; 
these abundances overlap the best-fitting regions from 
Madhusudhan et al. 2011) were also explored, under the 
assumption of chemical equilibrium, with similarly neg- 
ative results. Finally, we constructed a dayside model 
by dividing the planet's hemisphere into 10 concentric 
regions centered on the substellar point. These regions 
were modeled separately, each receiving flux along the 
line-of-sight to the star. The outgoing intensities dur- 
ing secondary eclipse, along the line-of-sight to the ob- 
server, were integrated to produce a dayside spectrum 
following Barman et al. (2005). The predicted limb-to- 
substellar temperature structure implies a horizontal as 
well as vertical temperature inversion (i.e., along a path 
from the substellar point to the terminator at a con- 
stant radius or pressure the temperature decreases, then 
increases again), but the surface- integrated fluxes from 
this model did not significantly differ from the default 
one-dimensional model. Our main conclusion from this 
modeling exercise is that it is very unlikely that strong 
narrow band emission could be reproduced by a mod- 
ification to the thermal profile alone, given the strong 
constraints placed on the SED by the broad band pho- 
tometry. 

As discussed above, the default, ID, model predicts a 
steep temperature inversion. A direct consequence of this 
inversion is a swath of narrow emission lines that form in 
the inverted atmospheric layers (grey-shaded spectrum 
in Fig. 6). While the narrow band filter does not encom- 
pass the complete CO band, it does cover many of the 
predicted narrow emission lines including the strong 3-1 
band. In the default model, these lines are very bright 
relative to the continuum but are too narrow to signifi- 
cantly impact the photometry. The model exercise above 
weakens the case for flux increases in the pseudocontin- 
uum, leaving the lines as a potential explanation for the 
deep narrow-band secondary eclipse. If the line fluxes 
are increased by a factor of 2 to 4, the narrow band pho- 
tometry can be reproduced. In this ad-hoc scenario, the 
lines would be extremely bright compared to those of 
the default model (see comparison in Fig. 6 inset). In 
the default atmosphere model, CO is the fourth most 
abundant molecule throughout most of the atmosphere 
including much the inversion region, as shown in Fig. 7. 
Previous studies have shown that local thermodynamic 
equilibrium (LTE) is achieved for CO in the atmospheres 
of isolated late-type stars (Ayres & Wiedemann 1989; 
Schweitzer et al. 2000) , meaning that the CO lines should 
map the gas temperature. However, Barman et al. (2002) 
concluded that Na was well out of LTE in the upper at- 
mosphere of the modestly irradiated planet HD 209458b, 
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resulting in strong emission cores for the Na D line pro- 
files. Nevertheless, the LTE assumption has not been 
tested for molecules in highly irradiated atmospheres. If 
such strong, narrow emission lines were present in the 
atmosphere of a hot Jupiter, they should be easily dis- 
cerned by future observations. Note that while CH4 also 
has a strong band head at this wavelength, (as observed 
in the telluric transmission spectrum and in spectra of 
L dwarfs; cf. Hinkle et al. 2003; Gushing et al. 2005), at 
the high temperature of WASP-12b atmospheric condi- 
tions would have to be much farther from equilibrium 
for significant CH4 opacity to be apparent. However, 
a careful non-LTE study is beyond the scope of this pa- 
per and further observational confirmation of the narrow- 
band emission is needed. 

The claim of strong, non-LTE CH4 emission in 
HD 189733b (Swain ct al. 2010), which was based on 
single-slit IRTF/SpeX spectroscopy, was disputed on the 
basis of contamination by telluric effects (Mandcll et al. 
2011). However, the types of spectroscopic contamina- 
tion discussed by Mandell et al. (2011) do not apply to 
our narrowband relative photometry. Changes in telluric 
absorption are common mode over our narrow field of 
view and should be removed when we divide the flux from 
WASP- 12 by the comparison star flux. Variable telluric 
emission should be removed by the combination of global 
sky frame subtraction and local background subtraction 
in the aperture photometry process (Section 2.2). As dis- 
cussed in Section 2.3, though the telluric absorption fea- 
ture intersecting our narrowband filter could explain the 
systematic photometric ramp (seen shortly after WASP- 
12 crosses the meridian; Figure 2 and Sec. 2.3), it seems 
unlikely that telluric variations could strongly corrupt 
the secondary eclipse depth while maintaining such a 
consistent eclipse duration and time of center (Sec. 3.3). 

Our analysis of the observations to date support a 
planet with little to no redistribution of flux to the night 
side, consistent with Cowan et al. (2012). We also find 
that photometric observations are well-reproduced by a 
black body and are not yet sufficiently precise to justify 
the use of more complicated models. If WASP-12b has 
a near-isothermal photosphere, then secondary eclipse 
data will be poorly suited to reveal significant compo- 
sitional information. Other highly irradiated giant plan- 
ets have also been observed to host nearly isothermal 
infrared photosphere, including the similarly hot WASP- 
18b (Nymcyer et al. 2011) and perhaps the cooler TrES- 
2b and TrES-3b (O'Donovan et al. 2009; Croll ct al. 2010; 
Cowan & Agol 2011). If WASP-12b does indeed largely 
radiate like a black body, previous conclusions about 
composition and the the presence or absence of a temper- 
ature inversion are significantly weakened. If our narrow 
band flux measurement is confirmed at higher precision 
and this flux is produced by emission lines, there may 
be hope for high-resolution spectroscopic studies to infer 
the planet's atmospheric composition and thermal pro- 
file (e.g. Barnes et al. 2007; MandeU ct al. 2011; Rodler 
et al. 2012; Brogi et al. 2012). 

7. CONCLUSIONS 

We have presented a deeper-than-expected secondary 
eclipse (0.45 % ± 0.06 %) of the very hot Jupiter WASP- 
12b in a narrow band centered at 2.315 fim. The 
planet's brightness temperature at this wavelength is 



3640 K±230 K, only marginally consistent with WASP- 
12b's equilibrium temperature of 2990 ± 110 K (Cowan 
ct al. 2012). Our precision is lower than expected because 
of an unanticipated systematic trend affecting both sky 
background and stellar photometry, but wc are able to 
exclude data affected by this trend. The duration and 
timing of the eclipse we measure from these data are 
consistent with a circular orbit and with previous mea- 
surements (Hebb ct al. 2009; Croll et al. 2011; Campo 
ct al. 2011; Cowan et al. 2012). 

Using NIK photometry and high-resolution spec- 
troscopy, we find that Bergfors-6, a previously identi- 
fied object only 1" from WASP-12 (Bergfors ct al. 2011, 
2012), is an M dwarf star with T^ff =3660t^^; K. If this 
object is an unresolved binary with two components of 
equal mass it could lie at the same distance from Earth as 
does WASP-12. However, Keck/NIRSPEC spectroscopy 
shows no evidence for binarity. If single, it likely lies 
closer to Earth than does WASP-12. Adaptive optics 
imaging on large-aperture telescopes will be necessary 
to conduct the proper motion studies necessary to dis- 
criminate between these two scenarios. If WASP-12 and 
Bergfors-6 are gravitationally bound, further simulations 
(e.g., Ibgui et al. 2011) should be undertaken to de- 
termine whether Kozai interactions with an object with 
Bcrgfors-6's characteristics could have caused WASP- 
12b's inward migration and, through tidal pumping, have 
inflated the planet's radius (Bodenhcimer et al. 2001). 

Bergfors-6 has heretofore passed unnoticed in previous 
transit and occultation analyses, which has caused the 
measured depths of these past events to be underesti- 
mated. We use our constraints on Bergfors-6 to infer 
and correct for the dilution of these past observations, in 
several cases increasing depths by > la. Thus WASP- 
12b is rather hotter and slightly larger (by 1 — 2%) than 
previously reported. These changes emphasize the im- 
portance of high-resolution imaging surveys in the vicin- 
ity of newly discovered transiting planets. 

The ensemble of dilution-corrected secondary eclipse 
measurements suggests that WASP-12b's atmosphere is 
largely isothermal across the pressures probed by eclipse 
observations (Figs. 6 and 7), with a photospheric tem- 
perature of roughly 3000 K. This result implies that pre- 
vious claims of a high carbon to oxygen ratio for this 
planet (Madhusudhan ct al. 2011) are not yet justified 
by the current photometric data. Further observations of 
the planet's 4.5 /xm secondary eclipse depth is certainly 
warranted to resolve the discrepancy between previous 
results at this wavelength (Campo et al. 2011; Cowan 
et al. 2012), as our modeling efforts indicate that achiev- 
ing such a dramatically lower temperature in the 4.5 /im 
channel would require a major unexpected change in 
the opacity source(s) across this bandpass. Regardless, 
our narrowband measurement alone excludes the mod- 
els used to infer this high C/0 ratio at > The lack 
of absorption at a wavelength where CO, a dominant 
species in any atmospheric model, should exhibit strong 
absorption is further evidence for a near-isothermal pho- 
tosphere. Thus secondary eclipse observations are ill- 
suited to determine WASP-12b's atmospheric composi- 
tion, and ultimately transmission spectroscopy may be a 
more successful approach in pursuit of this goal. 

WASP-12b is clearly an unusual object, and further 
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observations are clearly warranted. Aside from the need 
for additional IRAC 4.5 /im occultation photometry as 
described above, any or all of narrowband photometry 
(from the ground or, if available, using HST/NICMOS), 
single- or multi-object spectroscopy (Swain et al. 2010; 
Bean ct al. 2010; Mandeh et al. 2011; Berta et al. 2012), 
or perhaps by high-resolution phase curve spectroscopy 
(Barnes et al. 2007; Brogi et al. 2012; Rodler et al. 2012) 
could be of great utility. Finally, a better measurement of 
the planet's three-dimensional shape is also highly desir- 
able, especially given the apparent disagreement between 
the degree of prolateness inferred by Cowan ct al. (2012) 
and Swain et al. (2012). If WASP-12b is substantially 
prolate, three-dimensional models, ideally coupled with 
a general circulation model of the planet's atmospheric 
dynamics, may also provide important clues toward un- 
raveling the mystery of WASP-12b's atmospheric struc- 
ture and composition. 
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